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ABSTRACT 


The research tasks completed during this reporting period include an 
experimental study of the effect on the spread of flames of the turbulence 
intensity of an opposed air flow, and a theoretical analysis of the concurrent 
spread of flames over thin fuels. Both studies are in our opinion important 
contributions in the study of the fame spread process. The results of the 
experimental study show that the flame spread process is significantly affected 
by the flow turbulence intensity for flames spreading over both thin and thick 
fuels. For a fixed flow velocity, the spread rate decreases as the turbulent 
intensity is increased. This appears to be due to the turbulent convective 
cooling of the fuel surface and gas in the vicinity of the flame front. Also 
observed is that extinction of the flames occurs at lower velocities as the tur- 
bulent intensity increases. The results of the theoretical analysis, which are 
in good agreement with previous experimental measurements, give detailed infor- 
mation about the flame structure and mechanisms of flame spread. It is shown 
that the flame length and consequently the rate of flame spread, are strongly 


dependent on the interaction between the pyrol, .s and burn-out fronts. 
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1. INTRODUCTION 

The research performed during this reporting period forms part of an 
on-going research program aimed at the study of the mechanisms controlling the 
spread of fire in the presence of a forced flow of air, or other oxidizing gas. 
The oxidizing flow interacts with the fire and therefore affects its rate of 
spread. The final objective of the research program is to provide information 
about the spread of fire that can be used in the development of models of fire 
growth in rooms and of test methods of the flammability of materials. 

The tasks conducted during this period include experimental observations 
of the effect on the opposed flow flame spread of the flow turbulence intensity, 
and numerical modeling of the concurrent spread of flames over thin fuels. We 
feel that both works are important contributions in the study of the flame 
Spread process. 

In most practical cases the spread of fire occurs under turbulent flow con- 
dit.-1s. Flow turbulence can affect both the heat transfer and the chemical 
kinetic mechanisms that control the flame spread process, both by convective 
cooling and gas mixing. The determination of the role of these mechanisms is 
necessary to understand and predict the spread of fire under realistic con- 
ditions. With the exception of the works of Orloff, et al. (15th Int. Symp. 
Comb. p. 183, 1975) and Saito and co-workers (lst Int. Symp. Fire Safety, p. 

75, 1986; 24th Nat. Heat Trans. Conf. 1987) for the upward turbulent spread of 
flames, no systematic experiments have been performed to date to determine the 
effect of turbulence on the flame spread process. This is the objective of this 
work. The work consists of a parametric study where the velocity and turbulence 
intensity of an air flow are varied to study their effect on the rate of flame 


spread over thick and thin fuels. 


With regard to the theoretical task reported here, very few analyses have 
been published to date of the concurrent spread of flames over a thin com- 
bustible. The works of Markstein and de Ris (14th Int. Symp. Comb. p. 1085, 
1973) and Fernandez-Pello (Comb. Flame, 31, p. 135, 1978; Comb. Flame, 36, 

p. 63, 1979) provide only approximate analytical solutions of the problem, which 
can only describe partially the experimental observations. The present work 
provides a more accurate description of the process by removing some of the cri- 
tical assumptions made in previous analyses. The results describe and agree 
well with the eeenimen te observations. 

In the following section a summary is provided of the above studies. 

Papers prepared for publication describing these studies are included in this 


report in the form of appendices. 


II. RESEARCH PROGRESS 

II.1 Experiments of Flame Spread in an Opposed Turbulent Flow 

An experimental study has been conducted of the influence on the spread of 
flames over the surface of thick and thin solid combustibles of the turbulence 
intensity of an air flow opposing the direction of flame propagation. The 
experiments are conducted in a facility specifically designed to perform flame 
Spread experiments consisting of a small scale combustion tunnel and supporting 
instrumentation based primarily on optical and thermocouple measuring methods. 
The turbulence intensity in the tunnel test section is varied by means of grids 
and perforated plates of different sizes placed at the exit of the tunnel con- 
verging nozzle. The velocity and turbulence intensity of the flow field are 
measured with a one component Laser Doppler Velocimeter. The flame spread rate 
is obtained visually timing the pass of the flame front through marks placed at 


fixed distances along the combustible surface. 


II,1.1 Thick PMMA Sheets 

The measurements of the flame spread rate over PMMA sheets as a function of 
the opposed air flow velocity are shown in Fig. 1 for several values of the tur- 
bulence intensity. It is seen that the spread rate first increases sharply as 
the flow velocity is increased and then decreases. The initial increase appears 
to be the result of a transition from a mixed flow type flame spread to a forced 
flow one. At low flow velocities, normal buoyancy tends to lift the flame from 
the fuel surface, which results in reduced heat transfer from flame to fuel and 
paneeecently in a slower spread rate. The decrease of the flame spread rate 
with the opposed flow velocity at higher air velocities is the result of gas 
phase chemical kinetics effects. As the opposed flow velocity increases the 
residence time of the fuel vapor/air mixture ahead of the flame decreases. For 
constant chemical time this results in a slowdown of the combustion reaction at 
the flame leading edge that sustain the propagation of the flame. This effect 
is observed to occur at all turbulence intensity levels. 

The effect of turbulence intensity on the flame spread rate can also be 
deduced from the results of Fig. 1. This is presented in Fig. 2 for three value 
of the flow velocity. It is seen that for all flow velocities the flame spread 
rate first increases and then decreases as the turbulent intensity increases. 
The maximum occurs at approximately 6% of the turbulent intensity. This result 
may reflect the differences in the boundary layer that take place during the 
transition from laminar to turbulent flow. For PMMA, heat conduction through 
the solid is an important flame spread mechanism. Thus, the thinning of the 
boundary layer that occurs as the flow changes from laminar to turbulent should 


result in an increase of the flame spread rate because the heat transfer to the 
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solid increases. At low turbulence intensities this effect appears to be more 
important than the potential enhancement of the fuel convective cooling that 
Should occur as the flow becomes more turbulent. Another possible mechanism 
that should be considered is an enhancement by the flow turbulence of the fuel 
vapor and air mixing at the flame leading edge, which could result in a stronger 
gas phase reaction and consequently in a faster flame spreading velocity. 
Although we think that the former mechanism is the more likely to occur, there 


is not enough available information at this point to resolve the issue. 


11.1.2 Thin Paper Sheets 


The results of the measurements of the variation with the air flow velocity 
and turbulence intensity of the flame spread rate over thin paper sheets are 
presented in Fig. 3. In this case, the flames do not propagate (extinguish) in 
gas flows with velocities larger than 1.5 m/s. The low flow velocities at which - 
the tests had to be conducted limited the maximum range of turbulence inten- 
Sities attainable to 15%. From the results of Fig. 3 it is seen that for all 
turbulence intensities the spread rate decreases as the opposed gas velocity 
increases. The decrease in the low velocity region, however, is small. This 
again appears to be the influence of buoyancy on the flame spread process at low 
flow velocities. This influence is not as noticeable as in the case of PMMA 
probably because with paper the flames are smaller and consequently less 
buoyant. The stronger decrease of the flame spread rate with the flow velocity 
observed for the larger air velocities appear to be the combined effect of con- 
vective cooling and weakening of the chemical reaction (smaller residence time) 


that take place as the opposed flow velocity increases. 
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FIGURE 4 


An interesting result obtained in these experiments is the effect of flow 
turbulence intensity on the extinction of the flame. This is shown in Fig. 4 
where the gas velocity at which extinction is observed is presented as a func- 
tion of the turbulence intensity. It is seen that flame extinction occurs at 
lower flow velocities when the turbulence intensity is increased. This result, 
although it is in general agreement with the other flame spread results, is 
somewhat unexpected since turbulence is often viewed as a means to enhance fuel 
combustion processes. 

A paper reporting this result will be presented at the 1988 Fall 
Meeting/Western States Section of the Combustion Institute. It will also be 
submitted for public-tion in Combustion and Flame. The paper is included in 


this report as Appendix A. 


II1.2 Model of Concurrent Flame Spread Over a Thin Charring Solid 


A theoretical model has been developed of the spread of a flame over the 
surface of a thin combustible solid in a gas flow moving in the direction of 
flame propagation. The combustion reaction and the solid pyrolysis are modeled 
respectively with a one-step, second order, and with a zero order Arrhenius 
reaction. The charring of the solid is modeled by treating the solid as a two 
component material, the combustible component that gasifies, and the char that 
remains occupying the initial solid volume. To simplify the mathematical solu- 
tion of the problem, a Hagen-Poiseuille flow is prescribed, thus eliminating the 
need. for using the momentum conservation equations. The transient, reactive, 
gas phase balance equations of energy and species coupled at the interface to 
the solid phase energy equation are solved numerically to describe the evolu- 


tion of the flame and soot properties, and through them to predict the variation 


with time of the flame, pyrolysis and burnout distances and their dependence on 
the flow velocity. 

The model is applied to the description of the spread of flames over paper 
sheets in a concurrent air flow. In Fig. 5 the predicted variation with time of 
the flame, pyrolysis and burnout distances are presented for several flow velo- 
cities, The experimental data of Loh and Fernandez-Pello (lstl Int. Symp. Fire 
Saf. Sci., 65, 1986) are included for comparison. In Fig. 6 the predicted 
dependence on the flow velocity of the steady state pyrolysis spread rate is 
also presented with the comparison with the experiments. It is seen that the 
predictions of the model agree well with the experiments. Particularly impor- 
tant is the description of the transition of the flame spread process from an 
initial accelerative stage to a final steady state. The transition is deter- 
mined by the interaction between the progress of the pyrolysis and the burn-out 
fronts, and their effects in the flame length. According to boundary layer Aric 
ses, the model predicts a spread rate that increases, although weakly, with the 
flow velocity, while the experiments show a spread rate that first increases at 
low flow velocities and then becomes constant as the flow velocity is increased. 
Some reasons for this potential disagreement are discussed in the work. Finally 
the predicted flame extinction limit due to increased concurrent velocity is 
very close to the experimentally determined value. 

This work was performed in collaboration with Dr. Di Blasi of Naples 
University, Italy. A paper reporting the results has been accepted for publi- 
cation at the Proceedings of the 22nd Symposium International on Combustion, 
held in Seattle, Washington, August 1988, It is included in this report as 


Appendix B. 
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ABSTRACT 


An experimental study is presented of the influence on the spread of flames 
over the surface of a solid combustible of the turbulence intensity of a gas 
flow opposing the direction of flame propagation. In the study, the flow velo- 
city and its turbulent intensity are varied to observe the effect of these flow 
parameters on the rate of flame spread over thick PMMA and thin filter paper 
Sheets. Also observed is their effect on the extinction of the flame. The 
results show that for thick PMMA sheets the spread rate first increases and then 
decreases as the turbulent intensity is increased. The presence of a maximum 
appears to be the result of the flow transition from laminar to turbulent that 
influences the heat transfer process from flame to fuel. Through inter- 
ferometric photographs it is Shown that the decrease of the spread rate with the 
turbulence intensity is primarily due to the convective cooling of the region 
near the flame front. For thin paper sheets the spread rate always decreases as 
the turbulence increases. The observed differences between the thick PMMA and 
the thin paper results are probably due to their different mechanism of heat 
transfer ahead of the flame. It is also found that extinction of the flame 
occurs at lower flow velocities for higher turbulence intensities. This result, 
which is unexpected, suggests that convective cooling and/or flame stretch, and 
no turbulent mixing, are the major effects of the flow turbulence on the extinc- 


tion of a spreading flame. 


iP 


INTRODUCTION 


In most practical cases the spread of fire occurs under turbulent flow con- 
ditions. Flow turbulence can affect both the heat transfer and chemical kine- 
tics mechanisms that control the flame spread process. In opposed flow flame 
Spread, turbulence could produce an increased convective cooling of the unburnt 
fuel surface and gas region ahead of the flame. Also, it could increase the 
conductive heating of the solid due to the enhanced heat transfer in the burning 
region. The former effect should result in a decrease of the flame spread rate 
as the flow turbulence increases. The latter in an increase. Similarly, tur- 
bulence mixing could affect the flame spread process by either favoring the 
mixing of fuel vapor and oxidant and thus enhancing the combustion reaction, or 
by diluting the fuel vapor mixture with oxidant which could weaken the reaction. 
Determination of the role of these different mechanisms is necessary to 
understand and predict the turbulent spread of fire. 

With the exception of the work of Orloff et al, [1] for upward turbulent 
flame spread over large PMMA sheets, and of Saito and co-workers [2,3] for 
upward turbulent spread of flames over PMMA and wood sheets, no systematic 
experiments have been performed to date to determine the effect of turbulence on 
the flame spread process. This is the objective of this work. The work con- 
sists of a parametric study of the spread of flames over thick and thin fuels in 
a turbulent air flow opposing the direction of propagation. The parameters 
varied are the flow velocity and its turbulence intensity. The effect of these 


parameters on the rate of flame spread and on the flame extinction is reported. 
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Experimental Installation and Procedure 


The experimental apparatus consists of a small scale combustion tunnel, 
designed to conduct combustion experiments under varied flow conditions, and the 
supporting instrumentation, based primarily on optical and thermocouple 
measuring methods, A schematic diagram of the installation is shown in Fig. 1. 
The tunnel has a settling chamber 0.89 m long with a square cross section 0.3 m 
by 0.18 m, that is connected to the tunnel test section by a converging nozzle 
with a reduction area of 5.5 to 1 in 0.3 m. The tunnel test section is 0.6 m 
long and has a rectangular cross section 13 mm wide and 76 mm high. The side 
walls are made of 6 mm Pyrex glass, and the floor and ceiling sides are made of 
25.4 mm thick Marinite sheet. The tunnel is mounted horizontally on a three 
axis positioning table. The positioning of the table is accomplished with lead 
screws actuated with computer controlled stepping motors. The system permits 
positioning the test section with respect to the stationary optical instrumen- 
tation at 0.1 mm spacing. 

For the thick fuel experiments the specimen is mounted flush in the 
Marinite floor with its upstream edge placed 15 mm from the exit of the tunnel 
convergent nozzle. In the present work the fuel specimens are 12.7 mm thick 
PMMA (Rohm and Mass, Plexiglas G) sheets, 76 mm wide by 0.3 m long. The thin 
fuel specimens are 0.33 mm thick Whatman filter paper, 76 mm wide by 0.3 m long. 
The paper sheets are mounted in a metallic frame by inserting them in metallic 
spikes placed on the sides of the frame [4]. The paper is positioned in the 
middle of the test section to permit unrestricted burning over both surfaces. 
The paper sheets are dried in an oven for twenty-four hours prior to their use 
in the experiments, 

The gas flow in the tunnel is supplied either from centralized compressed 


air or from bottles of compressed oxygen and nitrogen. The individual flow 
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rates are controlled and measured with critical nozzles. Turbulence is 
generated by means of grids, or perforated plates, placed at the exit of the 
tunnel converging nozzle. A prescribed turbulence intensity is obtained through 
a combination of flow velocity and plate blockage ratio. As a consequence high 
turbulence intensity could not be obtained at very low flow velocities. The 
distribution of turbulence intensity through the test section is characterized 
by an initial region, approximately 40 mm where the turbulence decays sharply, 
followed by a 0.2 m region of practically constant turbulence, and a final 
region where the turbulence decays again due to entrainment effects at the test 
section exit. The gas velocity and turbulence intensity are measured with a one 
component LDV (TSI) operating in the forward scattered mode. The laser source 
is a Spectra Physics Ar.-Ion laser with a 2 W total output power. The signal is 
processed with a TSI Doppler processor and transferred to a micro-computer 
(IBM/AT). 

The experimental installation also includes a laser interferometric system 
that can be operated either as a Mach-Zehnder Interferometer or as a Holographic 
Interferometer. In this work the former arrangement was used to provide quali- 
tative information about the effect of the flow turbulence on the flame and 
thermal layer structure. In addition to the interferometer, the installation 
contains a thermocouple based set-up for gas and solid phase temperature 
measurements, 

The rate of flame spread is measured by timing visually the pass of the 
flame over previously ruled marks on the fuel surface. For PMMA sheets the 
flame spread process is initiated by igniting the top downstream edge of the 
sheet with a propane torch. The filter paper sheets are ignited with an 
electrically heated nichrome wire placed in close contact with the paper at its 


downstream edge. 
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RESULTS 


1. Thick PMMA sheets. 

The measurements of the flame spread rate over PMMA sheets as a function of 
the opposed air flow velocity are shown in Fig. 2 for several values of the tur- 
bulence intensity. Because of the difficulty of producing high turbulence 
intensity at low flow velocities, no data was taken at those conditions. The 
data correspond to averages of three different tests. The spread rate was 
measured in the test section region that has an approximately constant tur- 
bulence intensity. 

The results of Fig. 2 for low turbulence intensity show that the spread 
rate first increases sharply as the flow velocity is increased and then 
decreases. The initial increase appears to be the result of a transition from a 
mixed flow type flame spread to a forced flow one. At low flow velocities, nor- 
mal buoyancy tends to lift the flame from the fuel surface, which results in 
reduced heat transfer from flame to fuel and consequently in a slower spread 
rate. The effect of buoyancy on the flame and thermal layer location can be 
observed by comparing the interferograms shown in Fig. 3 for flow velocities of 
0.2 and 0.5 m/sec. It is seen that for 0.5 m/sec the thermal layer is thinner, 
and the flame slightly closer to the surface, than for 0.2 m/sec, which confirms 
the above argument. The decrease of the flame spread rate with the opposed flow 
velocity at hfaner air velocities is in agreement with the previous results of 
Fernandez-Pello et al. [5]. As inferred in that work, and later predicted 
mene catty’ [6], the decrease is the result of gas phase chemical kinetics 
effects. As the opposed flow velocity increases the residence time of the fuel 
vapor/air mixture ahead of the flame decreases. For constant chemical time this 


results in a slow down of the combustion reaction at the flame leading edge that 
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sustain the propagation of the flame. This effect is observed to occur at all 
turbulence intensity levels. At flow velocities larger than 2 m/s the flame 
spread rate is so slow that the PMMA surface recedes, producing a stagnation 
region that acts as a flame holder. Thus, for gas velocities larger than 2 m/s 
the flame propagation is really a mass burning process rather than a flame 
spread process and is not considered here. 

The effect of turbulence intensity on the flame spread rate can also be 
deduced from the results of Fig. 2. This is presented in Fig. 4 for three value 
of the flow velocity. It is seen that for all flow velocities the flame spread 
rate first increases and then decreases as the turbulent intensity increases. 
The maximum occurs at approximately 6% of the turbulence intensity. This result 
may reflect the differences in the boundary layer that take place during the 
transition from laminar to turbulent flow (Fig. 5). For PMMA, heat conduction 
through the solid is an important flame spread mechanism [7]. Thus, the 
thinning of the boundary layer that occurs as the flow changes from laminar to 
turbulent should result in an increase of the flame spread rate because the heat 
transfer to the solid increases. At low turbulence intensities this effect 
appears to be more important than the potential enhancement of the fuel convec- 
tive cooling that should occur as the flow becomes more turbulent. Another 
possible mechanism that should be considered is an enhancement by the flow tur- 
bulence of the fuel vapor and air mixing at the flame leading edge, which could 
result in a stronger gas phase reaction and consequently in a faster flame 
spreading velocity. Although we think that the former mechanism is the more 
likely to occur, there is not enough available information at this point to 
resolve the issue. 

The results of Figs. 2 and 4 also show that once in the turbulent flow 


regime, the flame spread rate decreases as the turbulent intensity increases. 
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This effect seems to be the result of convective cooling of the fuel surface and 
gas at the region in the vicinity of the flame leading edge. The sequence of 
interferograms of Fig. 5 shows this clearly. It is seen that as the turbulence 
intensity is increased the thermal layer fluctuates strongly, sometimes becoming 
very thin near the flame front. This indicates the inducement by the turbulent 
eddies of cold air from the free stream toward the fuel surface, which results 
in os a cooling of the fuel surface and of the gas phase chemical reaction. 
Both are effects that tend to decrease the rate of flame spread. 
2. Thin paper sheets 

The results of the measurements of the variation with the air flow velocity 
and turbulence intensity of the flame spread rate over thin paper sheets are 
presented in Fig. 6. In this case, the flames do not propagate (extinguish) in 
gas flows with velocities larger than 1.5 m/s. The low flow velocities at which 
the tests had to be conducted limited the miximum range of turbulence inten- 
Sities attainable to 15%. From the results of Fig. 6 it is seen that for all 
turbulence intensities the spread rate decreases as the opposed gas velocity 
increases. The decrease in the low velocity region, however, is small. This 
again appears to be the influence of buoyancy on the flame spread process at low 
flow velocities. This influence is not as noticeable as in the case of PMMA 
probably because with paper the flames are smaller and consequently less 
buoyant. The stronger decrease of the flame spread rate with the flow velocity 
observed for the larger air velocities is in agreement with the results of 
Fernandez-Pello et al. [5]. As it was discussed in that work, this appears to 
be the combined effect of convective cooling and weakening of the chemical reac- 
tion (smaller residence time) that take place as the opposed flow velocity 


increases. 
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The influence of the turbulence intensity on the flame spread rate is also 
deduced from the results of Fig. 6. It is seen that in this case the spread 
rate always decreases as the turbulence intensity increases, with no eee 
difference during the transition from laminar to turbulent as observed with 
PMMA. This is possibly the result of the different mechanisms of heat transfer 
that control the spread of the flame for PMMA or paper. For thin paper sheets, 
heat transfer through the gas ahead of the flame is the dominant mechanism of 
heat transfer during the spread of the flame [7]. Thus, the thinning of the 
boundary layer during the transition from laminar to turbulent flow should have 
a secondary effect in comparison with the enhancement of the convective cooling 
caused by the flow turbulence. It is therefore understandable that for flames 
spreading over paper sheets the simple appearance of flow turbulence could 
result in a decrease of the spread rate due to both convective cooling of the 
gas and cooling of the gas phase chemical reaction. 

An interesting result obtained in these experiments is the effect of flow 
turbulence intensity on the extinction of the flame. This is shown in Fig. 7 
where the gas velocity at which extinction is observed is presented as a func- 
tion of the turbulence intensity. It is seen that flame extinction occurs at 
lower flow velocities when the turbulence intensity is increased. This result, 
although it is in general agreement with the other flame spread results, is 
somewhat unexpected since turbulence is often viewed as a means to enhance fuel 
combustion processes. The fact that in this case the flow turbulence tends to 
facilitate the flame extinction indicates that convective cooling of the reac- 
tion zone, and/or flame stretching effects, are the dominant mechanisms in the 
flame extinction process, and that the enhancement of the fuel/air mixing by the 


turbulence has a secondary influence on the gas phase reaction. 
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CONCLUDING REMARKS 


The detailed parametric study conducted in the present work of the spread 
of flames in an opposed turbulent flow has shown that turbulence intensity has a 
Significant effect on the rate at which the flame spreads, and on the conditions 
leading to its extinction. The study has introduced a whole new aspect of the 
Flame spread problem that is of special interest given the fact that most prac- 
tical fires take place under turbulent flow conditions. The nature of the 
results indicate that the flow turbulence may affect both the heat transfer and 
the chemical kinetic mechanisms that control the spread of the flame. Although 
throughout this paper phenomenological arguments are given in an attempt to 
explain the different experimental observations, the data acquired is not suf- 
ficient to prove with certainty that the arguments are correct. 

In order to fully understand the effect of turbulence on the flame spread 
process it appears that it may be necessary to study its influence in each one 
of the flame spread controlling mechanisms, individually. This could be done by 
conducting experiments under conditions where the flame spread process is con- 
trolled by only heat transfer effects (high oxygen concentration and low gas 
velocity) or by only gas phase chemical kinetic effects (low oxygen concentration 
and/or high flow velocity). Also important would be the identification of the 
non-dimensional parameters that control the turbulent flow flame spread process, 
and the subsequent correlation of the results in terms of these parameters. 

Also interesting is the observation that an increase in the flow turbulence 
causes an earlier extinction of the spreading flame. The fact that this result 
is contrary to the more common observation that turbulence favors the combustion 


reaction by enhanced mixing suggests that the extinction of a spreading flame 


és 


oy 
may be controlled by mechanisms that are somewhat different than those in free 


flames. The elucidation of this result deserves further attention. 
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EEST#ORGE IGURES 


Schematic diagram of experimental installation. 


Flame spread rate over thick PMMA sheets as a function of the opposed air 
flow free stream velocity for several flow turbulence intensities. 


Interferometric photographs of flames spreading horizontally along PMMA 
sheets for opposed flow velocities of 0.2 m/sec and 0.5 m/sec. 


Variation of the flame spread rate over thick PMMA sheets with the flow 
turbulence intensity for several free stream velocities. 


Interferometric photographs of flames spreading horizontally along PMMA 
sheets in an opposed air flow of 1 m/sec and varied turbulent intensities. 


Flame spread rate over thin paper sheets as a function of the opposed air 
flow free stream velocity for several turbulence intensities. 


Dependence on the turbulence intensity of the opposed air flow velocity at 


which extinction of the flames spreading over thin paper sheets is 
observed. 
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MAIN FLOW VELOCITY (m/sec) 


FIGURE 2. 


FIGURE 3. 
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ABSTRACT 


A theoretical model is developed of tre sc-aaz° :* a “are cyer = : 
ofa thin combustible solia in a cas flew ~ovirg ‘n a: 
oagation. Many of the simplifying assumptions ~ade in crevious stuctes -3ve 
been removed. The combustion reaction and tre soiiz cytolysis are -saeied 
respectively with a one-step, second orde™, and witn 43 zero orter + 
rate. The charring of the solid is modeled Sy creating =7e solid 3s 2 two 
component material, the combustible componert ‘rat iasifies, anc tre crar tnaz 
remains occupying the initial solid volume. “9 simplify tne vatreratical 
solution of the problem a Hagen-Poiseuille ficw is orescrided, tmus eliminating 
the need for using the momentum conservation equatiors, The transient, c-eac- 
tive, gas phase balance equations of erexgy and scecies coupled at “re interface 
to the solid phase balance equations are solved sumecically to descrise tre ayo- 
lution of the flame and solid properties, 2nc =nroucn trem to credict the 
variation with time of the flame, cyrolysis ard Surn-out distances anda thei- 
dependence on the flow velocity. [re application of tne analysis to the spread 
of flames over paper sheets in a concurrent aim flow gives good agreement with 
previous experimental data. Particularly important is tne description of tne 
transition of the flame spread process from an initial accelerative stage to a 
final steady state. The transition is determined by the interaction between the 
progress of the pyrolysis and the burn-out fronts, and their effects on the 
flame length. According to boundary layer analyses, the model predicts a spread 
rate that increases, although weakly, with the flow velocity, while the experi- 
ments show a spread rate that first increases at low flow velocities and then 
becomes constant as the flow velocity is increased. Possible reasons for this 
disagreement are discussed in the work. Finally the predicted flame extinction 
limit due to increased concurrent velocity is very close to the experimentally 
determined value. 
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INT200UCT ied 

The concurrent, or flew assisted, mode of flame soread is of narticular 
importance in the development of a fire because of its -apid ard consequently 
hazardous character. [If tne combustible is thin, the oropagation of the 
"pyrolysis" or "burning" fronts is further characterized by the consumption of 
the combustible in the burning region upstream. This produces a propagating 
“burn-out” front. The interaction between the two fronts affects the charac- 
teristics of the flame downstream and consequently the flame spread rate. 

Very few analyses have been published of the concurrent spread of flames 
over a thin combustible. The works of Refs. [1-3] all provide approximate ana- 
lytical solutions of the problem by using simplifying assumptions. A detailed 
discussion of these assumptions is given in the review of Ref. [4]. The analy- 
ses provide explicit expressions for the rates of spread of the pyrolysis and 
burn-out fronts in the form of power law dependences on time that partially 
describe the initially acceleratory part of the spread process experimentally 
observed! >>, However, they cannot predict the transition to the final steady 
state of the process observed in the experiments, the phenomena related to the 
extinction of the flame or the effect of char formation. The studies of Ref. 
[6-8] do include char formation in the modeling of flame spread process, but 
their application to the present problem is limited because they are basically 
thick solid analyses that include many cestrictive assumptions. 

The present work provides a more accurate description of the problem by 
removing some of the critical assumptions made in the previous analyses of the 
problem. Finite Arrhenius type cates are introduced to describe the gas phase 


reaction and the pyrolysis of the combustible, thus semoving the flame sheet and 


She) 


ignition temperature aocroximations. Solid properties that vary from treir 35ri- 
ginal material to their cnar value are included to describe more accurately the 


combustible "burn-out" process. 


MODEL 

The model describes tre spread of a flame over the surface of a thin com- 
bustible solid in an oxidizer flow moving in the direction of flame propagation. 
A schematic of the model is shown in Fig. 1. A fuel sheet of specified length 
and thickness 1S placed in a wall of a channel where the gas flow is prescribed 
to be a constant properties Hagen-Poiseuille flow. The initiation of the flame 
spread process is achieved by igniting a small portion of the fuel sheet at its 
upstream edge by means of a pilot flame, which is active for a short period of 
time (.25 s). During this period the gas velocity is set equal to zero. Once 
ignition occurs, the pilot flame is turned off and the flow velocity set to its 
specified value, initiating the flame spread process. The magnitude of the 
pilot flame and ignition time are kept to a minimum to reduce the effects of the 
ignition process on the flame spread rate. As the combustible component of the 
solid gasifies the remaining char occupies the same volume of the initial solid. 

In the formulation of the mathematical model, several additional assump- 
tions are made to simplify the solution. In the gas the combustion reaction 15s 
modeled with a one-step irreversible chemical reaction of the form Ve + Vg * Vp, 
with a second order Arrhenius rate. With the flow field prescribed and constant 
thermal properties, the gas phase governing equations are: 


- energy balance 


ee dW el ret 
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In the solid phase analysis the solid is assumed to be thermally thin, and 
therefore the temperature and density distributions across its thickness are 
uniform, The pyrolysis of the combustible is modeled following the Anse ei 
proposed by Kung [9] for wood pyrolysis. The solid density and specific heat 
change during the pyrolysis orocess from their initial values to their final 
values for the char. The pyrolysis reaction is assumed to be described by an 
expression of the for S + Ve + ve» which implies a rate of material disap- 


pearance and product appearance given by: 


During pyrolysis the actual density of the material can be written in the form: 
Bret aout Ln 02 ¢/Deq) + Dee (4) 

The pyrolysis rate is assumed to be described by a zero order Arrhenius 

expression of the form: 


da/dt- = Aco, exp(-E./RT.) ) 


Combination of Eqs. (4) and (5) gives finally: 


dere dtass- even (oe g/Deg)] (eg - Oc) exP(-E./RT.) (6) 
With solid properties uniform across its thickness, the solid phase energy 
balance equation becomes: 
peCedTo/at = -do/dt (L + [- aeq/(0eq - ace la * 
ve ie 4 aes 
og/lo, Deploy oo Ae Ui aaa 43 gf ay! t + 


(3/93x) [a(k 7. )/ax] (7) 
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In £q, ‘7) tre cadiant: “lux: from the flame to the surface and from the surface 
to ambient have been neglected for simplicity vceasons, [¢ is acknowledged, 
however, that the error introduced by this simplification could 5e significant 
for relatively large scales and sooty flames. [t is worth noting that although 
the conductive term in the solid phase balance equation has been retained, it 
has been found that this term coes not play an important role in the propagation 
process and it can be dropped from tne equation without change in the results. 
This is to he expected since, for thermally thin solid fuels, the heat transfer 
to the unburnt material occurs mainly through the sas phase’?! 1 | 

The initial conditions for the above equations are that the combustible 
sheet and gas flow are at ambient conditions, i.e., at t = 0 the gas temperature 
and species mass fractions are: be “I osmieet LontaGonlo s Yog and Yr = Yrq and 
the solid temperature and density are T. = Tog and p. = ogg The boundary con- 
ditions are: iF * 1g: Yo * Yop: Ye = Yp = 0, Vy = Yip ate the: inflow, aaa 
outflow: aT ,/ax ai aY./ax 2 0 for i = I, P, 0, F; the upper wall of the chan- 
nel is inert and non-conducting: aT / ay = 0, aY;/ay = 0 for 1 | 2? eee 
the gas solid interface the boundary conditions that apply until burn-out are: 


Daya / yas im (Y;) for i = 1, 0, P, T. = toy wieweene 


°g g Caer 
pyrolysis mass flux given by m = Aco. exp(-E./R Wes In the fuel burned region 


Diateyey = Moe mal suo 


a zero normal gradient condition is applied to the temperature and species con- 

centrations, i.e., split boundary conditions are applied at the burn-out front. 
The system of non-linear, coupled differential equations (1-2) and (6-7) 

with their initial and boundary conditions are solved numerically. The finite 


difference method of Ref. [12] is used for the gas phase equations and a 
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splitting metnod, tnat combines a fourth order Punge-Kutta method and the rack. 
ward scheme, is used to solve the solid pnase ecuations, A variasta integracian 


grid is used in the x and y directions. 


RESULTS 

The above model is applied to describe the spread of a flame over tne sur- 
face of a thin paper sheet in a concurrent flow of air. ‘he paper sheet is 
37 m long and £33107 3m trick. nese dimensions were selected tne same as 
those used in the experiments of 2ef. [5] to compare the theoretical and ‘tne 


experimental results. The properties used for the paper aret3: 


i aah ol L.2o5=reuy kok }, Ket smleo25xl@ [kd sem Skye, SLs -n75380 kl /kgd, 
ogy = 046x107? [kg/m], o¢ = 0.03x10°> [kg/m], The data used for the pyrolysis 
reaction are A. = 1x10!9 cs74y, E. = 125,520 (kd/Kmole] which are in the range 
of values suggested in Ref. [14]. The gas properties are taken as those for air 
at ambient temperature. The one-step reaction model for the gas phase com- 
bustion refers to the monomer cellulose, and the kinetic and thermokinetic 
constants are: A, = 1x10!) [sm*/kg], &, = 112968 [kJ/kmole] (Ref. [13]) and at 
= 11857 [kJ/kg]. The value for the heat of combustion is taken equal to that 
reported in Ref. [15], and it gives good agreement for the predicted and 
experimental !® values of the gas phase temperature. The grid size in the x 
direction is very fine close to the ignition zone ‘ax = .12x1972m for a.length 
of 0.02 m), then it varies linearly from 0.12x107¢m to 0.94x1072m for a length 
of 0.05 m, and it becomes constant fax = 0.84x107°m) in the remaining part of 
the integration domain. For the y direction, a step equal to 0.5x1074m is used 
for a length of 0.12x107@m close to the fuel surface, a linear variation to 


0.55x1072 for y = 0.5x1972m and a constant distribution for the remaining part 


of the integration domain. The time integration step is equa! to axioms: Te 


43 


nas reen found that a furtner decrease in tre integration steps zoes ror chance 
- the accuracy level of the results. 

A characteristic example of the predicted flame structure and solid sroper- 
ties distributions is shown in Fig. 1. The calculations are for a flame 
Spreading over a paper sheet in a concurrent flow with a mean velocity of 9.57 
m/s at a time of 2 s, after the initiation of the flame spread process. In 
order to give a good representation of the flame structure details only a smal] 
part of the integration domain is plotted. The gas phase isotherms (Fig. la) 
and species equi-concentrations (Fig. 1b) distributions show that the flame, 
except at its upstream leading edge and tip, has a diffusional character with a 
relatively large reaction rate, Also, except at the flame leading edge, the 
normal temperature and species concentrations gradients are very large and the 
streamwise ones small. At the flame leading edge and tip regions, a quenching 
layer is formed near the fuel surface and there is penetration of fuel and oxi- 
dizer through the reaction zone. This indicates that in these regions the reac- 
tion rate is slow and that finite rate kinetics must be employed to describe 
them accurately. The solid temperature (Fig.la) and density (Fig. 1b) distribu- 
tions show the relationship between the gas and the solid phase processes. A 
significant increase in the solid temperature does not occur until the arrival 
of the flame tip. From then on the temperature increases steadily until the 
solid starts to pyrolyze at which point the temperature levels off marking the 
pyrolysis front. The initiation of fuel pyrolysis is Dest observed from the 
density profiles of Fig. 1b. The pyrolysis initiation is characterized by a 
decrease in the solid density, which continues until all the fuel is pyrolyzed 
and only the constant density char is left. The termination of the pyrolysis 


process is also characterized by a --arp increase in the solid temperature cue 


ag 


the accumulation of "eat sy the cnar. Either the constant char Jensity on snacs 
temperature increase determines the location of tre durr-out front, 

The evolution of the temperature field, as tne flame spreacs, is sresented 
in Fig. 2 for three consecutive times and the same flow conditions as above, 
The strong increase of flame length with time is the result of tne large amount 
of excess fuel vapor that 1s produced as the pyrolysis region, and the pyrolyzed 
fuel, increase. The ‘uel vapor is convected downstream from the pyrolysis front 
reacting with the air and consequently extending the flame in the downstream 
region, The results of Fig. 2 also show the propagation of the upstream leading 
edge of the flame. This propagation is linked to that of the burn-out front, 
and consequently to the solid phase processes. Under extinction conditions, 
however, the position of the leading edge would be controlled by gas phase kine- 
tic mechanisms. Under these conditions the burn-out front propagation rate 
could be faster than that of the pyrolysis front which could result in the 
extinction of the flame. 

The predicted evolution of the solid temperature and density profiles 
during the spread of the flame for the above stated conditions is shown in 
Fig. 3. These profiles provide information about the initiation of the solid 
heating process and pyrolysis, and of the depletion of the combustible com- 
ponent. This information is used to determine the location and progress of the 
pyrolysis and burn-out fronts and from them to calculate theirs rates of spread. 
The location of the pyrolysis front is determined from the temperature profiles 
by defining a reference or “pyrolysis front temperature", a method that is com- 
monly used in the experiments>’>. In this work this reference temperature is 
taken as 618 K as in the experiments of Ref. [5]. The location of the burn-out 


front is easier to determine because when all the combustible component of the 
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solid #6 gasified, its density cecomes that 3f tre char ard tre temneratiyre 
increases sharply as tnere is ro fuel gasification. 

Results as those presented in Fig. 3 are used to produce tne predictions 35f 
Fig. 4 for the variation with time of the pyrolysis and burn-out distances under 
several flow conditions. The experimental data reported in Ref. [5] are also 
included in the figure for comparison purposes. To reduce the effect of the 
fuel ignition process, the distances and time are given with reference to the 
peepectige values once ignition has been achieved. It iS seen that the model 
predicts well the experimental data with a spread process that is initially 
accelerative and that tends towards a steady State as the flame spreads over the 
solid surface. This is the result of the variation with time of the length of 
the pyrolysis region that initially increases very fast as the pyrolysis front 
propagates and the burn-out front still remains stationary. Once the burn-out 
front starts to propagate, the rate of increase of the pyrolysis length 
decreases as the burn-out front accelerates. When later the rates of spread of 
both fronts become of the same value, the pyrolysis length and consequently the 
rate of flame spread tends to a constant value. [In the same figure, the flame 
distances have been plotted. The flame tip is defined as the location of the 
downstream edge of the local fuel to oxygen equivalence ratio iso-line of value 
1. The predicted values for the flame length are also very close to the experi- 
mental ones. Better quantitative agreement between theory and experiments could 
be obtained for a given flow velocity by varying some of the parameters 
controlling the pyrolysis process such as the pre-exponential factor and activa- 


tion energy of the pyrolysis reaction, heat of evaporation, char density, etc. 
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It is felt, rowever, that the required computer time is not justifiable given 
the overall approximations of tre model and s--e of the uncertainties in tre 
measurements>. 

A difference between the model and the experimental “esults is the pre- 
dicted and measured dependence of the flame spread rate on the flow velocity. 
From the results of Fig. 4 it is seen that tne model predicts a spread rate that 
increases, although weakly, with the flow velocity, while the experiments Snow 3 
spread rate that first increases at low flow velocities and then becomes 
constant as the flow velocity is increased°, 

The predicted dependence on the flow velocity of the steady state spread 
rate of the pyrolysis front, deduced from the distances reported in Fig. 4 by 
differentiation with respect to time, is presented in Fig. 5. Also the experi- 
mental results of Ref. [5] have been included for comparison purposes. It is 
seen that the model predict a weak increase of the pyrolysis front spread rate 
with the flow velocity. A similar dependence can be deduced from Fig. 4 for the 
flame and burn-out fronts. Although these dependences can be explained with 


phenomenological descriptions? 


of the flame spread process, they disagree with 
the experimental observations of Ref. [5]. There are several reasons thay may 
explain these discrepancies. one is the type of flow prescribed in the present 
model which cannot predict the growth of the velocity boundary layer and con- 
sequently cannot predict accurately the surface heat flux distributions. 
However, the models using the boundary layer approximation for the flow field 
also predict similar dependences, which seems to indicate that this is not the 


primary problem. Another reason is the model simplified description of the 


filter paper pyrolysis and degradation process, which could result in inaccurate 
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predictions of tre Surn-out spread rate. The descciation of tre gas phase ~ea¢ 
tion is another potential source of error, altnough a sensitivity analysis of 
the gas phase reaction Arrhenius parameters effects on tre spread sate did not 
Show major differences in the results. Finally a fact that should be kept into 
account is that the conclusions described in Ref, £5] may not be totally 
accurate because of the experimental difficulties encountered and tne smal] 
range of flow velocities tested, 
As in the experiments, the extinction limit (fig. 5) iS characterized sy a 
discontinuous drop in the flame spread rate from a finite value to zero. In tre 
simulations it has been observed that a gas velocity of 2.18 m/s is the maximum 
value for which the spread of flames in air can occur. For higher velocities 
there is no propagation after ignition. The predicted limit, which is very 
close to that observed in the experiments (2 m/s), depends on the parameter 


values used to describe the combustion in the gas phase and the solid fuel 


pyrolysis. 


CONCLUDING REMARKS 

The present model provides a detailed description of the flame structure, 
as it spreads, and predictions of the flame, pyrolysis and burn-out distances 
variations that are in good agreement with the experimental observations. Also 
the determined extinction limit is very close to that experimentally observed. 
Particularly important is the description of the interact cntbereet the 
progress of the pyrolysis and burn-out fronts that results in the transition of 
the flame spread process from an initial accelerative stage to a final steady 
State. Another important contribution is the consideration of longitudinal dif- 


fusion during the initial stage of the flame spread and in the description of 
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the hurn-out process since tre soundary layer approximation cannct 4e 30piied as 
the upstream leading edge of the flame. The depencence of the ssread rate on 
the flow velocity is still a subject not fully resolved. Althoucn the experi- 
ments indicate that it becomes constant for relatively large flow velocities, 
the present and previous models of the problem predict a spread rate that 
increases weakly with the flow velocity. The dependence on the flow velocity is 
physically understandable since the spread of the pyrolysis front is controlled 
by mechanisms that are somewnat different than those of the burn-out front (mass 
burning rate), and it is the interaction between the spread of the two fronts 
that finally determines the flame spread rate. In addition there are uncertain- 
ties in the experimental data that suggest that perhaps the experimental conclu- 
sions are inaccurate and not the prediction of the model. This issue can only 
be resolved with more experimental information or more accurate models of the 
flame spread process. 

The results of this analysis also indicate that simplified versions of the 
present model could predict, it not in detail, certain aspects of the flame 
spread process. An example of such simplified models is the analysis of Saito 
et al’ for the upward turbulent spread of flames over charring materials. These 
simplified models would be particularly useful as input in models of fire deve- 


lopment in rooms, or in the application of tests to measure the flammability of 


materials. 
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NOMENCLATURE 


preexponential factor 

gas phase specific heat at constant pressure 
Solid phase specific heat 

gas phase diffusion coefficient 
activation energy 

heat of combustion ser unit mass fuel 
thermal conductivity 

length 

heat of pyrolysis 

pyrolysis mass flux 

molecular weight 

surface heat flux 

universal gas constant 

temperature 

time 

gas velocity 

mean gas velocity 

spread rate 

reaction rate 

coordinate parallel to the fuel surface. 
coordinate normal to the fuel surface 
mass fraction species 1 
stoichiometric coefficient 

density 

solid half thickness 


no) 
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SUBSCRIPTS 


combustible component in the solid 
burn-out 

char 

final char 

solid fuel that remains as char 
flame 

fuel 

gas phase 

solid fuel that vaporizes 

inert 

species 

Oxidant 

pyrolysis 

product 

solid phase 

ambient conditions 
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Fig. 1 - Schematic of model. Predicted flame structure and solid properties 
distributions for a mean flow velocity equal to 0.57 m/s, at t = 2 5 after isri- 
tion: a) gas isotherm field and solid temperature distributicn; 5) species 
Sete rercent cations fields and solid density distribution. The external] 
isotherm, for the gas phase, is 600 K ard tne increment is 300 K, The isolines 
POGenemsuelmassmt action correspond tony.4. 0.3, 00.2, 0.5) 0.901; ‘for the ox+- 


gen mass fraction the internal line is 9.05 and the increment is 0.05, 


Fig. 2 - Predicted evolution of the temperature field for a mean flow velocity 
Of 0.67 m/S, a).t = 1 5, 5) to= 1.75 us oc) tis.2.25 5, The external isotherm, 


for the gas phase, is 600 K and the increment is 300 K. 


Fig. 3 - Predicted evolution of the solid temperature and density during the 
spread of the flame in a flow with a mean velocity of 0.67 m/s from t = 45 to 


t = 15 s and with an increment of 1s. 


Fig. 4 - Predicted variation with time of the flame, pyrolysis and burn-out 
distances for several flow conditions. The experimental data of Ref. [5] are 


included for comparison. 


Experiments of Ref. [5]: 


gas velocity flame pyrolysis burn-out 


0.625 m/s ° o ¢ 
1.00 m/s X aN & 
2.00 m/s : CQ) Gi 

NOde]| | Sumner met LAMCya eel we PYCOlyS1 S408 OE) burn-out; 


mean gas velocity: a) 9.67 m/s, 6) 1 1/s, ¢) Pisa et Sd eany,S « 


og 


the Slow velocity of tne cysaMisis Scnaadl -are 
Figss5i< Predicted dependence on the y 


The experimental data of ef. [5] are included for comparison. 
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